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The gravitational waves and energy radiations from a spinning compact object with stellar mass 
in a circular orbit in the equatorial plane of a supermassive Kerr black hole are investigated in this 
paper. The effect how the spin acts on energy and angular moment fluxes is discussed in detail. The 
calculation results indicate that the spin of small body should be considered in waveform-template 
production for the upcoming gravitational wave detections. It is clear that when the direction of 
spin axes is the same as the orbitally angular momentum ("positive" spin), spin can decrease the 
energy fluxes which radiate to infinity. For antidirection spin ("negative"), the energy fluxes to 
infinity can be enlarged. And the relations between fluxes (both infinity and horizon) and spin look 
like quadratic functions. From frequency shift due to spin, we estimate the wave-phase accumulation 
during the inspiraling process of the particle. We find that the time of particle inspiral into the 
black hole is longer for positive spin and shorter for negative compared with the nonspinning particle. 
Especially, for extreme spin value, the energy radiation near the horizon of the extreme Kerr black 
hole is much more than that for the nonspinning one. And consequently, the maximum binging 
energy of the extreme spinning particle is much larger than that of the nonspinning particle. 
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I. INTRODUCTION 

It is a popular opinion that at the center of our Galaxy and many galaxies there are supermassive black holes JU-JU- 
And extreme mass ratio inspiral (EMRI): a compact stellar-mass object (black hole, neutron star or white dwarf) 
spiraling into the supermassive black hole in galactic nuclei is one of the most important sources of gravitational waves 
for NASA's Laser Interferometer Space Antenna (LISA) [3 ■ People estimated that LISA should see about 50 such 
events of EMRIs at distances out to redshifts z > 1 6] [7] [8] , and the signal-to- noise ratios of these sources are around 
10-100 [9]. Furthermore, the gravitational signals from EMRIs can also offer very useful information to study the 
mass distribution of galactic nuclei and the basic physics near the horizon of black holes, such as the no-hair theorem. 

For the motives mentioned above, we should have enough accurate and numerous waveform templates to match 
the upcoming observational data in the future. A very universal method to calculate the gravitational wave tem- 
plates is the black hole perturbation technology developed by Teukolsky [l(| [H| ■ The Teukolsky equation determines 
the (linearized) evolution of a perturbation to the Kerr spacetime. There are many published papers investigating 
gravitational radiations and orbital evolution. The Teukolsky equation is solved in these papers by post-Newtonian 
expansion ( 12][13]and references inside) and numerical simulation ([l4j]-[l8| and references inside). 

Usually, the study of EMRIs did not containe spins of the small mass particles. Because for the spinning particle 
in the Kerr spacetime, the dynamical equation is non-integrable, different from the spinless one. The motion of 
spinning particles is much more difficult to work out. Only a few papers considered gravitational waves from the 
spinning particles for some special orbits by using the Teukolsky equation. Mino, Shibata and Tanaka first investigated 
gravitational waves induced by a spinning particle falling into a rotating black hole [19]. Then Tanaka et al. researched 
the case of particle with small spin in circular orbits around a Kerr black hole[20|. 

On the other hand, when the spin is measured by terms of fJ,M, the spin parameter is 

where fi, M are the mass of the small object and central black hole respectively. People think that the spins of small 
objects are so tiny that they can be ignored without obvious influence. 

But small spin can produce considerable accumulated effect when we observe the secular evolution of EMRIs, and 
can also make the phase of gravitational waves different from the spinless case. The importance of spin was emphasized 
by Burko who considered a spinning compact object in a quasicircular, planar orbit around a Schwarzschild black hole 
[21] [22]. The effect of spin is at most marginally relevant for signal detection (23j . And we will see that the results of 



this paper can prove the spin of small body cannot be ignored. 

Furthermore, for the bigger value of spin, perhaps the first order perturbation method is not accurate enough, but 
the result can also offer some useful information on intermediate-mass-ratio inspirals (IMRIs). Especially, when spin 
achieves extreme value (S ~ 1), some very interesting phenomena appear. For example, in the Schwarzschild or Kerr 
background, researchers found chaotic behavior of extreme spinning particles [24| - J27| . Even without chaos, extreme 
spin can also produce some distinct orbital characteristic s 12711 . Obviously, gravitational signals of chaotic orbits would 
be much different from those of regular ones. References |28l] and (29j by using quadruple radiation formula calculated 
gravitational waveforms of chaotic systems. 

In the present paper, as the beginning of future research, we study the gravitation radiation of spinning particles 
in circular orbits on an equatorial plane by calculating the Teukolsky equation. In Ref. pol ]. the authors used the 
post-Newtonian expansion and assumed very small spin to keep the circular orbits. Differing from them, we calculate 
the Teukolsky equation numerically and make the particles be restricted on equatorial plane exactly. So we can get the 
results in a strong field where the small objects are near the black holes (r ~ M). The signals from such a strong- field 
area are easier for the detectors to catch. And also we can study gravitational radiation and orbital evolution when 
the value of spin becomes much bigger. 

As mentioned above, the strong spin of particle can produce chaotic orbits in some parameter range. But the 
circular orbits studied in this paper, can be determined completely and off course are nonchaotic. This ensures the 
stability of numerical results. In the future, we will discuss the gravitational waves and radiation reaction of chaotic 
systems in another paper. 

This paper is organized as follows. In the next section, we first use the Papapetrou equation to deduce the orbital 
parameters of the spinning particle in the equatorially circular orbit. Second, we introduce the Teukolsky-Sasaki- 
Nakamura formalisms briefly and our numerical method. In Sec. IV we present our results in detail. Finally, our 
conclusion and discussion are given in the last section. 

We use units G = c = 1 and the metric signature (—,+,+,+). Distance and time are measured by the central black 
hole mass M, and the energy of the particle is measured by fj,, the mass of itself. We measure the orbital angular 
momentum and spin of the particle in unit of fiM . 
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II. CIRCULAR ORBIT OF A SPINNING PARTICLE 

The equations of motion of a spinning test particle in a curved spacetime were given first by Papapetrou [30| , and 
then reformulated more clearly form by Dixon (3lj . In our units the can be written as 

dn^ 1 

^ = -(^*^P a S P + n a P a V, (3) 
— = -WK afh S a pP& + K a S a )v». (4) 

Where, is the four-velocity, p M is the linear momentum and the spin vector. The latter two must satisfy 

p v p u = -[i 2 , S V S V = S 2 , (5) 

where S represents spin magnitude. We define 



R-^vufi — n^/iv e pca/3j (6) 



and e paa P is the Levi-Civita tensor. 

Because of spin, the motion of particle does not follow the geodesic, then the p^ 1 is no longer parallel to . If 
following Dixon, choosing the rest frame of the particles center of mass, we can get one of the spin supplementary 
conditions, 

P% = 0, (7) 
and the relation between four-velocity and linear momentum 

where u v — p p '/ 1 p and spin tensor is 

S" v = ^ va(s S a u . (9) 

As with the spinless particle, there are two conservations, energy and total angular momentum which includes orbit 
and spin angular momentum 

e = -p t + \gt^S^, (io) 
Jz=P*-\g*wS>»'. (ii) 

Another constant, the Carter constant Q for the nonspinning particle, does not exit for a spinning one anymore. 

Usually, solving Eq. ([2HU) is difficult. By assuming the spinning particle moving in circular orbit on the equatorial 
plane of a Kerr black hole, the problem can be simplified greatly. In this assumption, the only nonzero component 
of spin must be s = S = —S/r. Then, using the condition of circular orbit dp r jdr = , combining Eq.Q and the 
relation (jHJ, after a long calculation, we get the angular frequency 



= d<t>/dt = =5 2 v — , (12) 

where s l presents the zth power of s. The upper sign refers to prograde and the lower to retrograde orbits and 

a = 2Mar, a± = -(3r 2 + 6a 2 )M, a 2 = {<iM 2 a + 3Mar), b Q = 4Mr 5 , (13) 

bi = -12Mar 4 , b 2 = 13M 2 r 4 , b 3 = -6M 2 ar 3 , b 4 = (9a 2 - 8Mr)M 2 r 2 , (14) 

c = 2{Ma 2 r - r 4 ), ci = -QMa(a 2 + r 2 ), c 2 = (4M 2 a 2 + 6Ma 2 r + 2Mr 3 ). (15) 
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This result coincides with the approximative expression Eq. (4.26) in Ref.© when s <C 1. Obviously, Eq. (|T2|) 
reduces to the angular frequency of the spinless particle if the spin s — 0, 



= d(f>/dt = ±- 



AfVa 



• 3 /2 ± oi/i/a ' 



(16) 



As mentioned before, because of the spin, the four-velocity is not parallel with linear momentum as before. From 
Eq. (JSJ), relation between them is 



v*=u* 



i ' = u ! 



M(3a 2 + r 2 )s 
" M(3a 2 + 2r 2 )s 2 



/j, 2 r 3 



3Ma(a 2 + r 2 ) 2 , 
23 3 u > 

3 Alas V 
/i 2 r 3 



(17) 
(18) 



We can clearly see that the difference between it M and w M is the s 2 order. And from the energy (I10|l and total angular 
momentum (fTTl) . we can get, 



Mt = 



£ (z-lE/h-XJx 



[l \ [1 + ZlS 

Jz , fz 2 J z /n + YE\ 



fi — z 2 s 



(19) 
(20) 



where, 



A/ 



^2 



M 



M 



(21) 
(22) 
(23) 
(24) 



Furthermore, the effective 4- velocity components 



■ .a 2 (r 2 +a 2 ) 2 
u =< + \\^ 



A 



^ + 3(1 



u*=ut+{ [i-^ + ^(i 



A 



(25) 
(26) 



where Mo' u o are the components without spin, and X, y just are the coefficients before s in Eq. (|19|20[) . 

For confirming a circular orbit, first the Kerr parameter a and spin magnitude S are fixed. Then we give a radius 
ro, from which 



dU(r) 



Or 



= 0, £ = t/(r ), 



(27) 



can determine the total angular momentum J z and energy E. Where U (r) is the so-called effective potential, 



D / D2 fi ,,2 y2 

U(J z ,a, S, r) = -J z + ^{—~ -)J| + 



(28) 
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and the related coefficients are 

„ rr, Ms 2 l9 , r .Mas.*,-, 

Z = VA(1 2/1- 2, 29 

A = [A + (r + a 2 M - z 2 s 2 - (1 (/x + z x s ) 2 Y 2 s 2 + ZFs, 30 

r r r 

B = (1 (u + zis 2 Fs — A + r 2 + a 2 (u - z 2 s) 2 Xs Z(l + XYs 2 ), (31 

r r r 

C = [A + (r 2 + a 2 )— - z 2S ) 2 XV - (1 - — + z lS f + —ZXs (32) 
r r r 

Finally the four-momentum components can be calculated from Eq. (|25l26[) . 

Now, we deduce the energy-momentum tensor of a spinning particle in curved spacetime. From Dixon's classical 
literature [32j, we can write down the energy-momentum tensor, 



T ali (x) = fi J dr 
= H I dr 



| * W ( ^ (T ' )) « (o (s, t)vV ( x , r) - V 7 [s^ a {x,T)v^{x,T) 



5^{x 


-z{t)) 




-9 


S^(x 


-z{t)) 




-9 


S (i \x 


-z(r)) 





fi I dr\ " ■^ r ^yj' u ^ _ \yn — » ^ s ^ (33) 

where z(t) is the world line of the spinning particle. For convenience, we write U a ^ = u^ a v^ + U^s. Considering 
there are only two independent components of the spin tensor in the circular orbit (29j. 

gtr _ _ S rt = ^ S r4, = _ g4>r = ^ (34) 

the nonzero components of are 

V%) = 1(4^ + T*^v* - T^utu 1 - T%u t u% 

Uf^-TttUtut + T^utut (35) 



I/ft = T^utu* + T^utuT - T r tt u^ - Y\ 



And the nonzero components of V ia " are 



V trt = A V** = -u^s; 

V rt * = ^(u t v t -u^)s. (36) 



Now, we project this energy-momentum tensor onto the Newman-Penrose null tetrad [32 



"•^,0,-^,, (37) 



" 2 VS S 

rh a — —= (ia sin 9, 0, £, — i(r 2 + a 2 ) sin#) . (38) 
v2 

It can be, 

T nn = n a n p T a0 , T nftl = n a m p T a ^ T™„ = m a m T a ^ (39) 
where p = (r — iacosO)^ 1 . So, 

T ab = liC ab (40) 
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where C a b can be expanded as 

l r A 2 nn a 2 A 2 sin 4 9 A aA 2 sin 2 »„ aAsin 2 # , % 

c nn = ^t 00 + t 11 + — - 2 — r 33 + 2-r 01 - 2 — - — r 03 - 2 — - — r 13 ], (4i) 

C nm = _^[ ^ Asing (r 00 + (o 2 + r2) sin 2 ^33) + msin0T Ol + AT «2_ 

2V2 E 

( r 2 + a 2 + a 2 S in 2 0)T 03 + £T 12 - * (r 2 + a 2 ) sin 0T 13 - a A sin 2 9T 23 } , (42) 

a-™™ = y [-a 2 sin 2 #T 00 + £ 2 T 22 - (r 2 + a 2 ) 2 sin 2 9T 33 + 2iaXsm9T 02 + 

2a{r 2 + a 2 ) sin 2 9T 03 - 2i£(r 2 + a 2 ) sin(9T 23 ]. (43) 
For circular orbits on the equatorial plane, the above equations can be greatly simplified. 

III. GRAVITATIONAL RADIATION AND RADIATION REACTION 

A. The Teukolsky-Sasaki-Nakamura formalism 

In this paper, we use the Teukolsky equation to calculate gravitational radiation and radiation reaction. The Teukol- 
sky formalism considers perturbation on the Weyl curvature scale ip^, which can be decomposed in the frequency 
domain [lfjj . 

/+oo 
dw^iW) - 2 SZ(dy m *e-^\ (44) 

-°° lm 

where p = — — ia cos 9) . The function Ri m u (t) satisfied the radial Teukolsky equation 

A2 i (i^r) - v{r)Ri - = (45) 

where ^ mw (r) is the source term, which will be given below, and the potential is 

K 2 + Ai(r - M)K n 
V(r) = ^ — + 8iu;r + \, (46) 

where K = (r 2 + a 2 )u ~ ma, A = Ei, n + a 2 uj 2 — 2amw — 2. The spin- weighted angular function _ 2 <S^"(#) obeys the 
following equation, 



1 d ( . A _ a Sj£ 



sm t 



smO d9 V d6 
The radial Teukolsky Eq. (H51) has the general solution 



. m 2 - Am cos 9 + 4 , 

(aw) cos 9 + Aaw cos 9 - 5 + E lm 

sin 9 



. 2 SZ = 0. (47) 



where the Rj% lLU ( r ) an< ^ Rfmu( r ) are two independent solutions of the homogeneous Teukolsky equation. They are 
chosen to be the purely ingoing wave at the horizon and purely outgoing wave at infinity respectively, 

R^ = B^A 2 e-^*, r^r + 

RL^B^ye^ + r-'BZ^re-^*, r -> oo; (49) 

and 

R%nu = DZie lp " + A 2 D^re"^* , r -> r + 

^ w =r 3 A^e---"*, r^oo, (50) 
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where k — lu — ma/2Mr+ and r* is the "tortoise coordinate". The solution (|48|) must be purely ingoing at horizon 
and purely outgoing at infinity. That is, 

iWr -> oo) = Z« w rV^*, (51) 
i?/ m .(r^r + ) = Z^A 2 e -^". (52) 

The complex amplitudes Z^^ are defined as 



J ' mU ^Bfl u J r+ A(r') 2 

rthole roo poo ( r f\ ( r f \ 

Zl -~ ^BZ^D-Jr Wf • (53) 

The particle in this paper is in circular orbit on the equatorial plane, thus the particle's motion is described only 
as the harmonic of the frequency Q^, (j 12[) . and defined 



io m = mO<£. (54) 



Then Z^°° are decomposed as 



Zlmw — Z h n S(uJ - UJ m ), 

Z% lu = Z% l 6(u>-u m ). (55) 
And the amplitudes Zi'°° fully determine the energy and angular momentum fluxes of gravitational radiations. 



B. Energy, angular momentum fluxes and waveforms 

From the Weyl scalar ip^, we can extract the two polarization components h + , h x of the transverse-traceless metric 
perturbation at r — > 00, 

^ 4 (r 00) -> ^(h+ —ihx.). (56) 



The energy flux in gravitational waves, from the Isaacson stress-energy tensor[39j is 

) , (57) 



KdAdtJ. 



^ 1 I fdh+V (dhy 



x 16tt \ V dt ) \ dt 
then we can get the fluxes at infinity as follow 

dEY ad ^\Z?* 



E 



The energy and angular momentum fluxes at the horizon can be calculated by the below formula ll], 

^V ad -v. 

-^J =2^ a '"> w • (59) 



Where the coefficient ai m is given as 



256(2Mr + )p(p 2 + 4e 2 )(p 2 + lGe 2 )^ 

= K^P ' (60) 
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with e = \JM 2 - a 2 /AMr+, and 

\Ci m \ 2 =[(A + 2) 2 + iamu m - 4a 2 w^](A 2 + 36maw m - SGa 2 ^) 

+ (2A + 3)(96a 2 w 2 „ - ASraauj m )lUuj 2 m {M 2 - a 2 ). (61) 

The total change in energy and angular momentum of the compact object is 

E = -(E H + E°°), j z = -(jj* + (62) 

From the energy loss and angular momentum loss of radiation, we can calculate the radiation reaction. There are 
two key points we can use to write down the orbit modification in an easy way: one is that circular orbits remain 
circular [35|-(37j; the other one is that S remains constant at first order 38] under adiabatic radiation reaction. The 
averaged rate of change of the radius is 

R (B 2 _C\j r> (B 2 _ Cy T 2 i ,.2(giy 

* = {k ^ %rffi '■>/«! r J- + - -We %,) 1 (63) 

It is obviously that spin can contribute to r in the first order. Finally, the gravitational waveform can be calculated 
by 

h+(0, 0, t)-ih x (9, & t) = £ ±-Z? m _ 2 S™- (6)e^-^ . (64) 



C. The Sasaki-Nakamura equation 

As mentioned in the above two subsections, for getting Zf^°°, we should integrate the homogenous version of Eq . (j4"5"j) . 
But there is a difficulty when one numerically integrates Eq. (|45|) due to the long-range nature of the potential V(r) in 
(14"5"1) . In order to solve this problem, Sasaki and Nakamura developed the Sasaki-Nakamura function X(r), governed 
by a short-ranged potential, to replace the Teukolsky function R(r) 34]. The Sasaki-Nakamura equation reads as 

F ( x )**™-U{r)X lmu = 0. (65) 



dr* 2 dr 

The functions F(r), U(r) can be found in Ref.[34j. The Sasaki-Nakamura equation also admits two asymptotic 
solutions, 

*bL=e-* r *, r^r + , (66) 
^l^C'/We-'+C/We-"*, r^oo; (67) 

and 

X^ u =C^ u P{ry^ + Ct u P{r)e-^\ r -> r + , (68) 
X^ = P{r)e-^\ r^oo, (69) 

where P(r), P{r) can be found in Refs.jioj. [ill ]. 

The solution of Eq. (1551 is transformed to the solution of the Teukolsky equation by 



r>H,oo _ 1 



A / yjr 2 + a 2 A dr ^r 2 + a 
The relations between the coefficients of the Sasaki-Nakamura function and the Teukolsky function are 



(70) 



A in 1 Aui 2 

oin Irrnj oholc D 00 (rj-i \ 



and the functions a, j3, r], and di muJ are given clearly in Ref.[34| or in Appendix B of |l z 
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D. The source term 



Now, we follow Ref. [3| to calculate the source term which is the left part of the Teukolsky equation (|45j) . For the 
circular orbit on the equatorial plane(# = 7r/2), 

= / dtA 2 {(A nn0 + A nrllQ + A r7irfl0 )6(r - r ) + d r [(A nrnl + A^ rTll )S(r - r )] + d^A^^r - r )]}. (72) 



All functions in the above equation are given as below 

-3^-1 



A nn 
AnfhO 
AfarhO 

Anrhl 
Afrirhl 



2pp p C nn rr trtc i n- t t el 
[L[L 2 S + 2iapL 2 S], 

2py/2p' 3 C nri 



SAt 

pSp^pC r - m 
St 



— - p- p) L 2 S + 1 1 — + p + p\ aS(p - p) 



iK 



A 



K 2 „ K „ 

£2 + 2*p-£ + ld r 



2/iV / 2/3 3 CVim r | 



SAt 
2pSp~3pCfj 

s7 

_pS£^p_ 
St 



[Lj5 + iap(p-p)5], 
iiT 



where 5 just is shorthand of _2<S£^[0(i)]. The projected energy-momentum tensor components are 

d 

Cab = U ab + iojV ab - imVt b + V r ab -^, 



(73) 
(74) 
(75) 

(76) 
(77) 

(78) 
(79) 



where ab represents nn, nfh, mm, and U ab , V 1 ^ can been calculated from U a ^ , V ia ^ using Eq. (|41|) - (j43|) . We can 
find that for the spinning particle there are three additional terms in the above equation compared to the nonspinning 
one. And 



L\S = auwa.es- ^ b k [(k - l)(fc + 2)] 1 / 2 _ i r fcm (0), 



k=l„ 



L[L\S= b k l( k -l)fc(fc + l)(fc + 2)] 1/2 Y km (6) + 2au sin 6 L\S - (aco sin 9) 2 S. 

fc=imin 



(80) 
(81) 



Where bi represents the coefficients of the eigenvector of £>, oYi m (8) and _iYj m ((9) are spin- weighted spherical har- 
monics. Of course, we should choose 8 = ir/2 in our case. 

Now, we can directly write down Z^, Zj^. Taking the source term (|72")l into Eq. ([5"3")) . then using Ea.(|5"5"]). we have 



iu m B" r 



4iujf n di mu] Bi- r 



5-^(ro), 



where. 



r-f/.OO, 



(ro) = i2^^(ro)[j4nrjO + 4m0 + ^mmo] 



,2 nff.oo 



dr 2 



4 - - 



mm2 • 



(82) 
(83) 

(84) 



E. Numerical method 



Now, we introduce simply our numerical method for calculating the gravitational radiation of the spinning particle 
inspiring Kerr black hole. First, we need determine the orbital parameters of the particle. Setting Kerr parameter a, 
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spin magnitude S and orbit radius, we solve Eg. (1271) to get energy and angular momentum of the spinning particle. 
Then based on the discussions in Sec. II, it is easy to get the angular frequency, four-velocity and linear momentum. 

Second, we numerically integrate the Sasaki-Nakamura equation to get Z^, Zj^ in the loop on I = 2, 3,4, • • ■ and 
to = — I, ■ ■ ■ ,1. But considering the symmetry of Z ;jr j°° about harmonics to, we only need a loop on to from 1 to I. 
Usually (based on the experience of many calculations done before) it is enough for obtaining satisfying accuracy to 
loop on 2 up to 8. 

Then from Eqs. ([58|) . (|59|) and (f62|) . we can calculate the energy and angular momentum loss due to the gravitational 
radiation. Finally, the orbit change of radiation reaction and waveform are calculated by Eqs. (j6"3")) and 

For checking the validation of our code, we compare our results with the post-Newtonian expansion at weak field 
without spin. They match very well. An an example, we give the comparison of the energy fluxes at infinity from 
a = to 1 while r = 25M in Fig. 1. In addition, we compare our results with the published data of Fujita and Tagoshi 
[lei when spin is zero, and find they have good agreement too. After this check, we think our numerical results are 
reliable. 



1.4x10" 

3* 1.35x10- 

w 8 1.3x10- 
•a 

1.25x10- 




0.0 0.2 0.4 0.6 0.8 

a/M 



FIG. 1. Comparison of the energy flux at infinity for orbit at r = 25M as a function of Kerr parameter a for I — 3, m = 1. 
Agreement between the numerical and post-Newtonian fluxes is good. 



IV. NUMERICAL RESULTS 

In this section, we present the numerical results of our calculations. We focus on how spin influences the gravitational 
waves. 



A. Orbital frequency with spin 

In Figs. ^ and ([3]), we show the shifts: Afi^ = f2^(5 ^ 0) — 17^(5 = 0) as a function of the orbital radius r while 
a = and 0.996 respectively. It is clear that when the direction of spin inverses, the orbital frequency is bigger than 
the nonspinning particle. This is because the total angular momentum is composed of the orbital angular momentum 
and the spin angular momentum of the particle. If spin has the same direction as the total angular momentum, the 
orbital angular momentum would decrease, whereas the orbital angular momentum must be added to counteract the 
spin. In addition, as the rotating angular momentum of the black hole increase, both the orbital angular momentum 
and the influence of spin on the angular momentum decrease. Usually, the absolute value of Afi^ increases as the 
orbital radius r decreases. But near the horizon of the Kerr black hole, perhaps there are some exceptions, for example, 
in Fig.©, the S = 10~ 2 one. 

And in Tables U and UH we list the detailed values of the frequency shifts at several orbital radiuses. It can be found 
that the influence of physical spin on the orbital frequency, as well as the frequency of gravitational waves, is small 
and can be negligible. But during the long evolution time of EMRIs ~ M 2 /fi, the accumulatively phasic error without 
considering the spin of the inspiraling object is considerable. The phasic shift due to spin should be calculated by an 
integration A</> = f^ v ° L \ AQ^(r(t))dt. This integration is difficult to calculate out, but we can do a simple estimate. 

For n/M ~ 10 -5 , the average frequency shift is about 10~ 7 (see Table Hi]) of the extreme Kerr black hole. During 
highly relativistic regime, inspiraling about ~ M/fi cycles, and the typical frequency ~ 10 -2 , the dephasing of 
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TABLE I. The orbitally angular frequency shift due to the spin for a — 0. 



s 


-lCT d 




-10~ 4 




-10~ 5 








10~ 4 




io-- 1 




T 


in 

— 1U 


1.50008 x 10" 


c> 


1.50001 x 10" 


7 


1.50000 x 10" 


8 


-1.50000 x 10" 


-X 


-1.49999 x 10" 


7 


-1.49992 x 10" 


(i 


r 


= 8 


2.92991 x 10" 


6 


2.92971 x 10" 


7 


2.92969 x 10" 


-8 


-2.92969 x 10" 


-8 


-2.92966 x 10" 


7 


-2.92946 x 10" 


6 


r 


= 6 


6.94527 x 10" 


6 


6.94453 x 10" 


7 


6.94445 x 10" 


-8 


-6.94444 x 10" 


-8 


-6.94436 x 10" 


7 


-6.94362 x 10" 


-6 


r 


= 4 


2.34426 x 10" 


5 


2.34380 x 10" 


6 


2.34376 x 10" 


-7 


-2.34374 x 10" 


-7 


-2.34370 x 10" 


6 


-2.34324 x 10" 


5 


r 


= 2 


1.87616 x 10" 


4 


1.87512 x 10" 


5 


1.87501 x 10" 


-6 


-1.87499 x 10" 


-6 


-1.87488 x 10" 


5 


-1.87384 x 10" 


-4 






TABLE II. 


The orbitally 


angular frequency shift ASl^ due to the spin for a - 


= 0.996. 




s 


-10 _J 




— 10 -4 




-lO"" 5 




10"" b 




10" 4 




io- 3 




r 


= 10 


9.65822 x 10" 


-7 


9.65768 x 10" 


-8 


9.65762 x 10" 


-a 


-9.65761 x 10" 


-y 


-9.65756 x 10" 


-8 


-9.65702 x 10" 


-7 


r 


= 8 


1.74151 x 10" 


-6 


1.74137 x 10" 


-7 


1.74136 x 10" 


-8 


-1.74135 x 10" 


-8 


-1.74134 x 10" 


-7 


-1.74120 x 10" 


-6 


r 


= 6 


3.61475 x 10" 


-6 


3.61431 x 10" 


-7 


3.61426 x 10" 


-8 


-3.61425 x 10" 


-8 


-3.61421 x 10" 


-7 


-3.61377 x 10" 


-6 


r 


= 4 


9.30694 x 10" 


-6 


9.30480 x 10" 


-7 


9.30459 x 10" 


-8 


-9.30454 x 10 


-8 


-9.30433 x 10" 


-7 


-9.30219 x 10" 


-6 


r 


= 2 


3.03510 x 10 


-5 


3.03302 x 10" 


-6 


3.03281 x 10" 


-7 


-3.03276 x 10" 


-7 


-3.03255 x 10" 


-6 


-3.03047 x 10" 


-5 


r = 


= 1.2 


2.55531 x 10" 


-5 


2.55116 x 10" 


-6 


2.55075 x 10" 


-7 


-2.55065 x 10" 


-7 


-2.55024 x 10" 


-6 


-2.54610 x 10" 


-5 



spinning and nonspinning particles is 



Acj) 



A M 2?r 
A^ — — 



2rr. 



(85) 



We can find that though the wave-frequency shift of spin is very small, during the evolution period in a highly 
relativistic regime, it can produce considerable (about a 1 cycle) phase error. From the above estimate, we think the 
role of spin in the EMRIs should be calculated to obtain the template of gravitational waveform accurately. 



0.002 



0.001 



9 0.000 



-0.001 



-0.002 




FIG. 2. The shift of the orbital frequencies between spin and nonspinning particles as a function of radius r while a = 0. From 



the top down, the spin magnitude S is —10 



10" 



10" 



10 respectively. 



B. Energy fluxes vs spin 



Now we show the results of energy fluxes of the spinning particles. In Figs. (QJ [5] and [5]), we list the energy fluxes of a 
particle orbiting the black hole with a — 0, 0.996 and —0.996 respectively at infinity and the horizon. All the orbital 
radiuses are r = 10. We can find the energy fluxes to infinity and the horizon being a function of the spin S looks 
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like quadratic functions. Usually, the dE^/dt decreases as the spin magnitude increases. This result coincides with 
the known relation between the energy fluxes and black hole rotation. In other words, the spin of particle together 
with the black hole's rotating angular momentum, if they are "positive" (the same direction with orbital angular 
momentum), will reduce the energy flux down to infinity, and if "negative", will add dEao/dt. At the same time the 
fluxes down to the horizon are no longer monotonic varying with spin; this is exhibited clearly in Figs HHHl except for 
the dE\ l (2, l)/dt for a = 0. This phenomenon is also analogous with the case of the dE^/dt vs the black hole's angular 
momentum a. An interesting phenomenon is that when a < 0, the function of dE^/dt of S is a concave function, but 
a convex function when a > 0. 

In FigJTJ we give function images of the gravitational wave luminosity vs the spin S, and we produce enough 
accuracy of the luminosity up through 1 = 8. All of the results we figured out should be enough to clear the effect of 
spin acting on the energy fluxes. The angular momentum flux, from Eqs. (|58l I59[) . is just the energy flux over fi^. It 
is not necessary to list the results of the angular momentum flux here. 



C. Energy loss vs spin 



Because of the gravitational radiation, the energy of a particle decreases , and so does the orbital radius. The 
total energy loss can be calculated by Eq. (|62|) . Through computing the energy- loss rates for different spin values at 
a certain radius, we can compare how the spin influences the energy loss of a particle due to gravitational radiation. 
For the convenience of our comparison, we introduce relative differences of the energy-loss rates, 

Aifo) = E(S ± 0) - E(S = 0) 
E E{S = 0) 

In Table Hill the values of AE(S)/E are listed in detail. 

It is assumed that LISA will observe gravitational waves of EMRIs at the typical frequency ~ 10~ 2 Hz and the total 
wave cycle is about N ~ 10 5 for one year. Thus, the relative difference of energy luminosity due to spin AE(S)/E 
required to establish the accuracy for the cycle AN < 1 must be < 10~ 5 in circular orbit cases [42|]. For more 
complicated orbits, such as zoom-whirl orbits, the requirement for accuracy would be stronger than for this circular 
one. From Table IIII1 we can clearly find for typical S = 10~ 5 that the relative differences of luminosity reach even 
exceed the above requirement. The results here argue again that the spin of a small object should be computed for 
the data analysis of LISA et al., and is consistent with the previous calculation of frequency shift in Sec. IV A. 

If we want to calculate the total energy radiation during the whole inspiraling process from a big r to the horizon, we 
need to integrate E step by step. But this calculation requires too much computer time, for example, for nonspinning 
particles, perhaps a CPU year or more. For spinning particles on circular or equatorial orbits, it would not add much 
cost [we need to calculate additional three terms in Eq. (|79)> ]. But for general orbits, it needs to numerically work 
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FIG. 5. The energy fluxes of I = 2 modes at infinity and the horizon as a function of the spin magnitude S. Where the Kerr 
parameter a = 0.996 and the orbital radius r = 10M. 
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x 1 0~ 7 _ x 1 




S/uM S/uM 

FIG. 6. The energy fluxes of I = 2 modes at infinity and the horizon as a function of the spin magnitude S. Where the Kerr 
parameter a = —0.996 and the orbital radius r = 10M. 




S/uM S/uM 

FIG. 7. The gravitational wave luminosity at infinity and the horizon (up through / = 8 to produce the accuracy) as a function 
of the spin magnitude S. Where the Kerr parameter a = 0.996 and the orbital radius r = 10M. 



out the Papapetrou Eq. (|2]|4j). By previous experiences [ISlH: we know that it adds about half cpu year compared 
to nonspinning cases. Fortunately, thanks to the theory which circular remains circular under adiabatic radiation 
reaction, the energy loss of the particle is just the energy difference of the particle at the circular orbit radius r and 
near the horizon. In Table ITVl we list the particle energy at radius r — 10M and r = 1.2M (very near the horizon), 
and give the energy difference between the two cases. We can find that the particle with positive spin will radiate more 
energy during the inspiral to the black hole than nonspinning one. At the same time, the particle with negative spin 
radiates less energy. Remember that the gravitational wave luminosity decreases when spin is positive and increases 
when spin is negative(see Fig-((7j)). This means the inspiraling time of positively spinning particle should be longer, 
and shorter for negative one. 

Especially, for the extreme spinning particle (S ~ 1), the radiated energy is much more than nonspinning one (see 
the first row in Table HV|) . This phenomenon only occurs when o « 1 and S is extreme (and the same direction as 
the rotation of black hole). It is known that the maximum fractional binding energy per unit rest mass is 42% for 



TABLE III. The relative differences of energy loss AE/E vs varied spin values for a — 0.996. 



s 


10"" 3 




io- 4 








-10~ 5 




-10~ 4 




-10~" 3 




r 


= 10 


-4.31993 x 10" 


4 


-4.32215 x 10" 


5 


-4.32303 x 10" 


(i 


4.32140 x 10" 


6 


4.32255 x 10" 


b 


4.32475 x 10" 


4 


r 


= 8 


-5.98264 x 10" 


4 


-5.98605 x 10" 


-5 


-5.98697 x 10" 


fi 


5.98834 x 10" 


fi 


5.98760 x 10" 


5 


5.99122 x 10" 


4 


r 


= 6 


-9.31789 x 10" 


4 


-9.32626 x 10" 


-5 


-9.32635 x 10" 


6 


9.32736 x 10" 


-6 


9.32800 x 10" 


5 


9.33652 x 10" 


4 


r 


= 4 


-1.84199 x 10" 


-3 


-1.84453 x 10" 


-4 


-1.84480 x 10" 


5 


1.84482 x 10" 


-5 


1.84510 x 10" 


4 


1.84764 x 10" 


3 


r 


= 2 


-9.76634 x 10" 


-3 


-9.60818 x 10" 


-4 


-9.59232 x 10" 


5 


9.58882 x 10" 


-5 


9.57298 x 10" 


4 


9.41434 x 10" 


3 


r = 


= 1.5 


-7.57872 x 10" 


-2 


-7.47114 x 10" 


-3 


-7.46037 x 10" 


4 


7.45798 x 10" 


4 


7.44723 x 10" 


3 


7.33930 x 10" 


2 



TABLE IV. The particle energy at r = 10M, r = 1.2M and energy loss vs varied spin values for a — 0.996. 



S 



1.0 
0.1 

10"" 3 
10 -5 


-io- 5 

-10" 3 
-0.1 



E(r = 10M) E(r = 1.2M) 



AE 



0.9504110109413 0.2966990902654 0.6537119206759 

0.9517720321309 0.6831909453624 0.2685810867685 

0.9519178894668 0.7337036132776 0.2182142761892 

0.9519193435776 0.7342411704276 0.2176781731500 

0.9519193582782 0.7342466042358 0.2176727540423 

0.9519193729484 0.7342520379715 0.2176673349769 

0.9519208270598 0.7347903703284 0.2171304567313 

0.9520657813190 0.7932472876811 0.1588184936379 



a = 1 and about 30% for a — 0.998 [43]. But we point out for extreme spinning particle, the binding energy is much 
larger than the nonspinning one. For example, for a = 0.996 and S = 1, the maximum fractional binding energy is 
about 77%. But we should take care that the physical value of spin is <C 1. Thus, this big binding energy here only 
has theoretical meaning. 



V. CONCLUSION AND DISCUSSION 



In this paper, we have calculated the gravitational waves emitted by a spinning particle in circular orbits on 
equatorial plane around a Kerr black hole. We solved the completely orbital parameters of the spinning particles: 
orbital frequency, energy and total angular momentum and nonzero four-momenta. We also give the relation between 
the four-velocity and linear momentum. 

Using the Teukolsky formalism of the Kerr black hole perturbation, we analyze the effect of spin on gravitation 
waves in detail. Our calculations were done in the strong-field area and can apply to extreme spin value. Perhaps 
this is the first time that the numerical results indicate the spin of a small body cannot be ignored for the future 
gravitational waveform data analysis . 

The main effects of spin presented in this paper are: a positive spin can decrease wave frequency and energy 
luminosity, but this result is reversed when spin is negative; the relations of energy flux of spin are close to quadratic 
functions; in the whole inspiral process, the evolution time of positively spinning particle is longer and the total energy 
radiated is larger; the overall phase difference compared to the nonspinning case achieves the critical value, 1 cycle; 
the maximum binding energy of extreme spinning particles is much more than for a nonspinning one. 

In future, it will be interesting to discuss the spin effect for more complicated orbits such as zoom-whirl orbits with 
a big eccentricity even with inclination. 
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